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ABSTRACT 

SRC variables are M supergiants, precursors to Type II supernovae, that vary in brightness with mod- 
erately regular periods of order 100-1000 days. Although identified as pulsating stars that obey their own 
period-luminosity relation, few have been examined in enough detail to follow the temperature and spectral 
changes that they undergo during their long cycles. The present study examines such changes for several 
SRC variables revealed by CCD spectra obtained at the Dominion Astrophysical Observatory (DAO) during 
2005-2009, as well as by archival spectra from the DAO (and elsewhere) for some stars from the 1960s to 
1980s, and Cambridge radial velocity spectrometer measures for Betelgeuse. Described here is our classifica- 
tion procedure and information on the spectral type and radial velocity changes in three of the stars. The 
results provide insights into the pulsation mechanism in M supergiants. 

Subject headings: Surveys; stars: individual: SRC variables: a Ori, a Her, S Per; stars: late-type. 



1. Introduction 

Type C semiregulars, or SRC variables, constitute 
a portion of the "forgotten" family of late-type vari- 
able stars located in the M supergiant region of the 
Hertzsprung-Russell diagram. The more irregular LC 
variables are close cousins, and both types are poorly 
studied, mainly because of their lengthy primary cycles 
of 200 — 900 days, although some also have superposed 
longer periods measuring thousands of days in duration. 
They are evolved stars, most being M supergiants, with 
a few bright giants and a few late K-type stars thrown 
into the mix. They have masses of order 15 — 20M Q , ef- 
fective temperatures below 4,000 K, and radii of at least 
several hundred Rq . Many are surrounded by large cir- 
cumstellar dust shells created by their ongoing mass loss. 
SRC variables are also believed to be the precursors to 
most Type II supernovae. 

Recent summaries of the properties of SRC variables 
were presented by Turner (2006), Rohanizadegan et 
al. (2006), and Turner et al. (2006), the last two studies 
in connection with the luminous M3 la variable BC Cyg. 
The stars appear to obey their own period-luminosity re- 
lation (Turner et al. 2006). Spectroscopic studies of a 
large sample of M supergiants were made by Levesque 
et al. (2005), including many of the SRC variables, and 
a study of period changes in BC Cyg reveals erratic 
changes of period in excess of what is typical of stars 
of late spectral type (Turner et al. 2009). It was demon- 
strated by Stothers (1969) forty years ago that pulsation 



can account for the primary periods of variability in such 
stars. But observational confirmation has been lacking, 
and detailed observational studies of the stars are few 
in number, mainly restricted to a limited spectroscopic 
and radial velocity survey by Joy (1942), the photomet- 
ric survey of BC Cyg by Rohanizadegan et al. (2006) 
and Turner et al. (2006), and a brief survey of a Ori by 
Gray (2008). 

In an attempt to remedy that situation, a spectro- 
scopic survey of northern hemisphere SRC variables was 
initiated in 2005 using the 1.85m Plaskett telescope 
of the Dominion Astrophysical Observatory (DAO), 
Herzberg Institute of Astrophysics. Examination of 
archival photographic spectra for bright SRC variables, 
some from the David Dunlap Observatory (DDO), and 
scans of many of the original DAO plates have also 
been made, as well as a search for unpublished spectro- 
scopic observations of the stars. Included are previously- 
unpublished radial velocity measures for Betelgeuse ob- 
tained by Griffin with the Cambridge radial velocity 
spectrometer. A comprehensive survey of the variables 
holds the promise of learning more about the basic prop- 
erties of an intriguing group of massive stars, objects 
that may explode as supernovae within a few hundred 
years. This study examines changes in spectral type 
and radial velocity in three SRC variables: a Orionis, 
a Herculis, and S Persei. It is part of a larger project 
designed to improve our understanding of the pulsa- 
tion mechanism in M supergiants and the evolutionary 
changes in SRC variables. 
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2. Data Reduction 

CCD spectra from the DAO were reduced with IRAF, 
first by applying bias corrections, trimming the spectra, 
and removing bad pixels with the CCDPROC package. 
Iron-argon arc spectra were used for wavelength calibra- 
tion, and flat fielding was done with the APFLATTEN 
package. The DO SLIT package was used to extract one- 
dimensional (1-D) spectra from the two-dimensional (2- 
D) CCD output, as well as to wavelength-calibrate and 
dispersion-correct the former. The spectra were then 
normalized using IRAF's CONTINUUM package, which 
fits a polynomial to the data to find the continuum, and 
then outputs the ratio of the input spectra to the fitting 
function. 

3. Spectral Classification 

Luminous M stars have been notoriously difficult to 
classify, primarily because most of the available spectral 
standards are SRC variables (Keenan & McNeil 1976). 
A fairly straightforward scheme of temperature classifi- 
cation was set up by Keenan (see Gray & Corbally 2009), 
based upon the visibility of specific band heads of TiO 
in blue-green spectra of the stars: A4954 at MO, A4761 
at Ml, A4804 at M2, A4584 at M3, A4626 at M4, A4352 
at M5, A4395 at M6, A4310 at M7, and distinct TiO ab- 
sorption longward of A4100 at M8. That same scheme 
was employed here, along with the use of a VaO fea- 
ture at A4389 that is a useful luminosity discriminant 
for warmer M supergiants. Table [1] summarizes the lines 
used to determine spectral type and luminosity class for 
this study. 



Table 1: Lines used for spectral classification. 



A (A) 


Species 


Behavior 


4861 


H/3 


Appears to strengthen with decreasing 
temperature when compared with 
nearby continuum, which is suppressed 
by the TiO band head. 


4953 


TiO 

band 

head 


Strengthens with decreasing 
temperature. A4953 and other band 
heads are the primary determinants of 
spectral type for M stars. The A4953 
band head is easiest to use. 


4389 


VaO 
blended 
with 
Fe I 


Strengthens with increasing 
luminosity when compared with 
nearby Fe I line at A4383 A. 



3.1. Radial Velocity Determination 

The IRAF package RVIDLINESwas used to calculate 
radial velocities from the CCD and photographic spec- 
tra. RVIDLINES uses the differences between observed 
wavelengths in a spectrum and rest wavelengths of the 
same lines to calculate a radial velocity. The user creates 



a list of rest wavelengths, then marks several lines in the 
spectrum and enters the rest wavelengths of the marked 
lines. The software then uses the user-created list to 
identify as many more lines in the spectrum as possible, 
and computes a velocity from the average wavelength 
shift, along with an associated uncertainty. When the 
software is run in heliocentric mode, it uses the observa- 
tory location keyword in the header along with the date 
and time information to apply a heliocentric correction 
to the computed velocity automatically. 

4. Results 
4.1. a Orionis 

Table 2: a Orionis observations. 



HJD 


Phase 


Vr fkm/s) 


Sp. Tvne 


Source 


2433937 


0.4036 




Ml. 5-2 lab 


DDO 


2437957 


0.9526 




M2 lab 


DDO 


2439921 


0.6167 




Ml. 5-2 lab 


DDO 


2441043 


0.2829 


+19.9 ±0.4 


Ml la 


DAO pg 


2441132 


0.4943 


+3.7 ±0.1 


Ml la 


DAO pg 


2441233 


0.7342 


+10.3 + 0.2 


M2 lab 


DAO pg 


2442458 


0.6439 


-21.9 + 0.8 


M3 lb 


DAO pg 


2442819 


0.5014 


-22.6 + 0.1 


M3 lb 


DAO pg 


2454041 


0.1562 


+23.5 




Griffin 


2454057 


0.1942 


+23.8 




Griffin 


2454079 


0.2464 


+23.3 




Griffin 


2454124 


0.3530 


+22.6 




Griffin 


2454186 


0.5023 


+18.5 




Griffin 


2454360 


0.9140 


+18.7 




Griffin 


2454391 


0.9888 


+22.3 




Griffin 


2454414 


0.0421 


+24.7 




Griffin 


2454421 


0.0587 


+24.9 




Griffin 


2454429 


0.0797 


+24.7 




Griffin 


2454443 


0.1110 


+24.7 




Griffin 


2454472 


0.1797 


+22.2 




Griffin 


2454490 


0.2246 


+21.8 




Griffin 


2454508 


0.2673 


+21.4 




Griffin 


2454523 


0.3029 


+21.0 




Griffin 


2454757 


0.8571 


+22.8 




Griffin 


2454779 


0.9091 


+25.1 




Griffin 


2454804 


0.9683 


+25.5 




Griffin 


2454845 


0.0678 


+26.6 




Griffin 


2454876 


0.1415 


+26.9 




Griffin 


2454897 


0.1911 


+26.3 




Griffin 


2454929 


0.2670 


+25.4 




Griffin 


2455114 


0.7050 


+25.1 




Griffin 


2455128 


0.7383 


+25.2 




Griffin 


2455167 


0.8307 


+23.9 




Griffin 


2455187 


0.8780 


+23.4 




Griffin 


2455214 


0.9443 


+23.2 




Griffin 


2455227 


0.9729 


+23.3 




Griffin 


2455258 


0.0487 


+25.2 




Griffin 


2455295 


0.1364 


+24.6 




Griffin 



A comparison of AAVSO photometric information on 
the star with newly-derived and available spectral type 
and radial velocity data indicates that a Ori reaches its 
latest spectral type and lowest luminosity near light min- 
imum. Radial velocity maximum (greatest photospheric 
recession) is reached near light maximum and the star 
reaches mean radial velocity (indicative of smallest di- 
mensions) near light phase 0.27. Table [2] summarizes 
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Fig. 1. — a Ori data. From top: visual light curve from 
AAVSO measurements, recent spectra, spectral types and 
temperatures as functions of phase, and radial velocity as a 
function of phase. 
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Fig. 2. — q Her data. From top: visual light curve from 
AAVSO measurements, recent spectra, spectral types and 
temperatures as functions of phase, and radial velocity as a 
function of phase. 



the spectroscopic results and Fig. [T] displays the avail- 
able data for a Ori. 



4.2. a Herculis 

A comparison of AAVSO photometric information on 
a Her with newly-derived and available spectral type 
and radial velocity data indicates that it reaches small- 
est dimensions at phase 0.12 following light maximum. 
There is scatter in its spectral classifications, yet it ap- 
pears to reach greatest luminosity at earliest spectral 
type. The star belongs to an optical triple. Table[3]sum- 
marizes the spectroscopic results and Fig. [5] displays the 
available data for a Her. 



4.3. S Persei 

A comparison of AAVSO photometric information on 
S Per with newly-derived spectral type and radial veloc- 
ity data indicates that it also reaches smallest dimen- 
sions at phase 0.11 following light maximum, and ap- 
pears to reach greatest luminosity and earliest spectral 
type near light maximum. Table 2] summarizes the spec- 
troscopic results and Fig. [3] displays the available data 
for S Per. 

5. Conclusions 

All three survey stars appear to reach their small- 
est dimensions following light maximum at photomet- 
ric phases 0.11-0.27. The photospheric temperature 
changes follow the changes in overall size, with a Ori 
and S Per reaching their coolest photospheric temper- 
atures near light minimum and greatest atmospheric 
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Fig. 3. — S Per data. From top: visual light curve from 
AAVSO measurements, recent spectra, spectral types and 
temperatures as functions of phase, and radial velocity as a 
function of phase. 



extent. Like S Per, a Ori reaches its highest photo- 
spheric temperatures just prior to reaching its smallest 
overall dimensions. Scatter in the spectral classifica- 
tions for a Her make it difficult to determine exactly 
when it reaches hottest photospheric temperature, but 
in all three stars greatest luminosity appears to be co- 
incident with such phases, implying that the T^ s term 
dominates the luminosity equation. Similar analyses are 
currently being performed on ~ 20 additional SRC vari- 
ables, which will permit us to determine if other mem- 
bers of the class display the same general characteristics. 



Table 3: a 


Herculis observations. 




HJD 


P llclSC 


V fkm/s") 




Source 


9497Q7S 


0.3787 




M5 lab-lb 


DDO 


2427984 


0.3905 




M5 lab-lb 


DDO 


2439293 


0.6529 


-65.9 ±0.7 


M4 la 


DAO pg 


2440033 


0.1089 




M5 lb 


DDO 


2441047 


0.1056 


-58.7 ±0.4 


M4 la 


DAO pg 


2442513 


0.9915 


-20.5 ±0.3 


M4 la 


DAO pg 


2454027 


0.6568 


-37.1 ±8.1 


M5 lab 


DAO CCD 


2455086 


0.7415 


-55.2 ± 19.8 


M6 lb 


DAO CCD 


Table 4: S Persei observations. 


HJD 


Phase 


V r (km/s) 


Sp. Type 


Source 


2453648 


0.7871 


66.6 ± 14.8 


M4 la 


DAO CCD 


2454379 


0.6929 


23.7 ± 13.5 


M5 lab 


DAO CCD 


2454749 


0.1514 


12.4 ± 15.9 


M3 Ia+ 


DAO CCD 


2455087 


0.5703 


20.7 ± 19.7 


M4.5 la 


DAOCCD 



(http://pollux.graal.univ- montp2.fr), operated at GRAAL 
(Universite Montpellier II - CNRS, France) with the support 
of the PNPS and INSU. 
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